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Abstract 



We revise the 7r"'"7r~ and qq annihilation mechanisms of dilepton production during relativis- 
' tic heavy-ion colUsions. We focus on the modifications caused by the specific features of in- 

CN . tramedium pion states rather than by medium modification of the /o-meson spectral density. 

The main ingredient emerging in our approach is a form-factor of the multi-pion (multi-quark) 
^ ' system. Replacing the usual delta-function the form-factor plays the role of distribution which, 

in some sense, "connects" the 4- momenta of the annihilating and outgoing particles. The dif- 
ference between the c.m.s. velocities attributed to annihilating and outgoing particles is a 
\ particular consequence of this replacement and results in the appearance of a new factor in the 

■ formula for the dilepton production rate. We obtained that the form-factor of the multi-pion 

. (multi-quark) system causes broadening of the rate which is most pronounced for small invariant 

\ masses, in particular, we obtain a growth of the rate for the invariant mass below two masses 

Q_i' of the annihilating particles. 

1^ , Since leptons do not interact practically with the highly excited nuclear matter, they leave the reac- 
' tion zone after their creation in a relativistic nucleus-nucleus collision without further rescattering. 
That is why, the dileptons {e^e~ and fi~^fi~ pairs) observed in high-energy heavy ion collisions se- 
^ • cure an excellent opportunity for obtaining information on the initial state and the evolution of the 
system created in a collision. The enhancement with an invariant mass of 200 -r- 800 MeV observed 
by the CERES collaboration Pl Ej in the production of dileptons has received recently a considerable 
attention and has been studied in the framework of various theoretical models (for the review, see 
Ref. jS]). It was found that a large part of the observed enhancement is due to the medium effects 
(see Refs. jHE] and references therein). Meanwhile, pion annihilation is the main source of dileptons 
which come from the hadron matter [HI d- That is why, the proper analysis of the dilepton spectra 
obtained experimentally gives important observables which probe the pion dynamics in the dense 
nuclear matter that exists at the early stage of the collision. The purpose of the present letter is 
to look once more on the tt+tt" annihilation mechanism of dilepton production from the hadron 
plasma by accounting the medium-induced modifications of the dilepton spectrum. In order to do 
this, we concentrate on the modifications which are due rather to intramedium pion states, than on 
the discussion of a modification of the p-meson spectral density. In accordance with our suggestions, 
the main features of a pion wave function follow from the fact that pions live a finite time in the 
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system where they can take part in the annihilation reaction. As a consequence, the off-shell effects 
give an appreciable contribution to the features of the annihilation process specifically in the region 
of low invariant masses. Moreover, if pions are the entities of a local subsystem, then the spatial 
structure of the pion states is far from a plane-wave one and this also gives the essential contribution 
to the features of the dilepton spectrum. 

To carry out the outlined program, we assume that the pion liquid formed after the equilibration 
exists in a finite volume, and the confinement of pions to this volume is a direct consequence of the 
presence of the dense hadron environment which prevents the escape of pions during some mean 
lifetime r. The same can be assumed concerning a hot system of quarks which are confined to 
a quark-gluon droplet. We sketch the proper geometry in Fig. 1. The small circle of radius R 
represents a subsystem of pions which is in a local thermal equilibrium and moves with the collective 
velocity v. If the fireball size is of order of i?o oc 6 -h 10 fm, then the mean size of the small 
subsystem is of order of R oc 1^5 fm. So, we assume the system of pions (quarks) produced in 
high-energy heavy-ion collisions is effectively bounded in a finite volume. The pion (quark) wave 
functions (px{x), where A is a quantum number, satisfy the proper boundary conditions and belong 
to the complete set of functions. For instance, the stationary wave functions may be taken as the 



solutions of the Klein-Gordon equation (V + k ) (fxi'^) 



0, where k"^ = E"^ — m?, which satisfy 
the Dirichlet boundary condition on the surface 5*: ip\{'x)\s = 0. For the box boundary, we get 
V9k(x) = \/^JV Y[^=i(^{Li — Xi)9{xi) sin (kiXi), where V = L1L2L3 is the box volume, A = k = 
{ki, k2, ks), and components of the quasi-momentum run through the discrete set ki = urii/ Li with 
rii = 1,2,3,... . For the spherical geometry, the normalized solutions are written as (fkimi'^) = 



9{R- 



■r) (2/r)^^^ Ji^i/2{kr)Yim{'&, (p) / RJi+3/2ikR) , where A = {k,l,m). Next, the field operators (p{x) 
corresponding to the pion field should be expanded in terms of these eigenfunctions, i.e. 

d^k 



0{x) 



(27r)32c<;k 



[a{k)Mx) + b'-{kM{x)] 



(1) 



where a(k) and 6(k) are the annihilation operators of positive and negative pions, respectively. 

On the other hand, the states corresponding to confined particles 
can be written in a common way as v^k(x) = a/ p(x) /V $k(x), where 




P 



1/2/ 



Fig. 1. Sketch of an expanding fire- 
ball. The small circle of radius R 
represents the subsystem of pions 
which is in a local thermodynamic 
equilibrium and moves with the 
collective velocity v. 



= nti [^(Li - Xi)e{xi)] for a box and p^/\x) = e{R - r) 
for a sphere, respectively. The function p(x) presents the infor- 
mation about the geometry of a reaction region or cuts out the 
volume where the pions (quarks) can annihilate. Hence, for the 
evaluation of S-matrix elements wave functions (pk{x) should be 
taken as the pion m-states once annihilating pions belong to fi- 
nite system. The amplitude of the pion-pion annihilation to a 
lepton pair in the first non-vanishing approximation is calculated 
via the chain tt+tt" ^ p — 7* — //, where the p-meson ap- 
pears as an intermediate state in accordance with the vector me- 
son dominance. The matrix element of the reaction is (out|S'|in) = 
- / d^xi d''x2 (p+, p_ I T [H]{x,) Hj{x2)] I ki, ka), where H[{x) = 
— e j^(x) y4'^(x) and H\{x) = — e j^(x) A'^(x). It is remarkable that 
the pion density p(x) appears as a factor of the pion current. Indeed, 



[x] 



V 



-i^{x) df, $" 



X 



(2) 



where the field operator <l>(a;) is defined in the same way as that in (Q) with the functions v^k(x) 
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replaced by $k(x). Because of this factorization, after the integration over the vertex x the density 
p(x) automatically cuts out the volume, where the tt+tt^ annihilation reaction is running. At the 
same time, this means that the density p(x) determines the volume of quantum coherence, i.e. just the 
particles from this spatial domain are capable to annihilate one with another and make contribution 
to the amplitude of the reaction. To obtain the overall rate, it is necessary then to sum up the rates 
from every coherent domain of the fireball. 

For the sake of simplicity, we assume that the pion states can be approximately represented as 
V^k(a^) = a/ p{x) /Ve~^^'^ . Here, p{x) is the 4-density of pions in the volume V (marked by radius R 
in Fig. 1), where pions are in a local thermodynamic equilibrium. In essence, this approximation 
considers just one mode of the wave function $k(x). We expect that this ansatz reflects the qualitative 
features of the pion states in a real hadron plasma, and, below, the main consequences of the 
contraction of particle states in the dense environment are considered. 

A simple calculation immediately shows that the ^-matrix element is proportional to the Fourier- 
transformed pion density p{x), i.e. (out 15*1 in) oc p{ki + k2 — p+ — p^), where ki and ^2 are the 
4-quasi-momenta of the initial pion states and p+ and p_ are the 4-momenta of the outgoing leptons. 
This means that the form-factor of the pion source p{k) stands here in place of the delta function 
which appears in the standard calculations, i.e. {2'n')'^6'^{K — P) ^ p{K — P), where K = ki + k2 and 
P = p+ + ]9_ are the total (quasi-) momenta of pion and lepton pairs, respectively. An immediate 
consequence of this is a breaking down of the energy-momentum conservation in the s-channel of the 
reaction, which means that the total momentum K of the pion pair is no longer equal exactly to the 
total momentum P of the lepton pair. The physical interpretation of this fact is rather obvious: the 
effect of the hadron environment on the pion subsystem which prevents the escape of pions from the 
fireball can be regarded during the time span r as the influence of an external nonstationary field. 
The latter, as known, breaks down the energy-momentum conservation. From now, the squared 
form-factor \p{K — P)p of the pion system plays the role of a distribution which in some sense 
"connects" in s-channel the annihilating and outgoing particles instead of 5-function. Indeed, the 
number N^^^ of produced lepton pairs from a finite pion system related to an element of the dilepton 
momentum space, reads 

(^^)^ .),^(^), (3, 

where 

/ dN \ f d^ki d^k2 „a,, , „ , ^ , „ , ^ , 

Here, Ei = ^/m?~+~^, i = 1, 2 for pions and i = for leptons, respectively. By the broken 

brackets, we denote the thermal averaging over the pion quasi-momentum space with the help of the 
thermal distribution function /th(-E). To obtain eq.Q, we represent the amplitude Ap = (outj^^^^ |in) 
of the reaction under consideration as 

Ap{ki,k2;p+,p-) = p{ki + k2-p+ - p^)Ao{ki,k2;p+,P-) . (5) 

We note that not only the form-factor p{K — P) contains information about the pion system. The 
amplitude Aq carries new important features as well, which are related to the violation of the energy- 
momentum conservation in the s-channel. Indeed, the pion-pion c.m.s. moves with the velocity 
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wk = K/ii'o; whereas the lepton-pair c.m.s. moves with the velocity vp = P/Pq- Hence, these 
two center-of-mass systems are "disconnected" now, that is why any quantity should be Lorentz- 
transformed when transferred from one c.m.s. to another. A particular consequence of this is the 



1 + i 



1 ( {P-KY 



in the formula for the dilepton production 



appearance of the correction factor 
rate (for details see [H]). 

Concerning the physical meaning of eq. Q, we note that one can regard it as the averaging of 
the random quantity ( d-^Kd^p ) '^i^h the help of the distribution function \p{K — P)\'^ centered around 
the mean value P. In this sense, the hadron medium holding pions in a local spatial region for 
some time, which is expressed as the local pion distribution p{x), plays the role of an environment 
randomizing the pion source. This randomization is a purely quantum one in contrast to the thermal 
randomization of the multi-pion system which is already included to the quantity { ^iRd-^p )- 

In order to transform the distribution of the number of created lepton pairs over the dilepton 
momentum space to the distribution over invariant masses, one has to perform additional integration 

using {dN^P^/d'T) from ©, i.e. (f^) = I^il^), where Pq = VM^ + P2 
together we obtain 



Taking this 



a 



3(27r) 



Ami 



1/2 



1 + 



2m2 



1 - 



M2 

4^2 X 3/2 



-f3Ko 



1 + 



- 1 



eiKo)eiK' - Ami) , (6) 



where F.,r{M^) is the p-meson form-factor, and we take the Boltzmann distribution /th(P) = exp(— /JE") 
with (3 = 1/T as inverse temperature. By the presence of the ^-functions (last two factors on the 
r.h.s. of ©), we would like to stress that the invariant mass of a pion pair = \/K^ is not less 

than two pion masses. On the other hand, possible finite values of the distribution ^^pr-^ below the 
two-pion mass threshold can occur just due to the presence of the pion system form-factor p{K — P). 

The factor in the square brackets on the r.h.s. of © is a correc- 
tion which is due to the Lorentz transformation of the quantity 
(ki — k2)^ from the dilepton c.m.s. to the pion-pion c.m.s. This 
factor gives a remarkable contribution to the dilepton spectrum for 
invariant masses below the two-pion mass value. Its influence is 
especially pronounced for e~^e~ production (see 0). 

To clarify the effects under investigation as much as possible, 
for particular evaluations we take as a model of pion system the 
Gaussian distribution of the particles in space and the Gaussian 
decay of the system of pions: 




p{x) 



exp 



2r2 2P2 



(7) 



Meanwhile, it can be another choice of the model function for the 
pion source. Indeed, one can choose for instance a geometry with 
sharp boundaries which are determined by the ^-functions. To 
elucidate this point we compare two formfactors (normalized to 
the unit volume) which correspond to the Gaussian distribution 

Pg{r) = e 2jj2 and to the 6'-function distribution pe(r) = ^(-R~|r|) (see Fig. 2). Only a slight difference 



Fig. 2. Comparison of the Gaus- 
sian and 0— function form-factors 

|p(q)IV^- 



4 



between these formfactors is seen and, therefore, the choice of pion source distribution doesn't affect 
much the dilepton production rate. In what follows, we dwell on the Gaussian distribution. 

To be closer to the quantities measured in experiment, we evaluate the rate related to the proper 
rapidity window 

2ttM— I dy / dP± 

J P I min 



where tanhy 
?(p) 



dM dy 

p3/p0^ p2 __ 



^ ymm 

+ (P2)2. 



d^x d'^P 



The results of evaluation of the production rates 
(i-Rg+''g_ / dMdy and dR^^l^_ / dMdy for electron-positron and muon-muon pairs, respectively, in pion- 
pion annihilation are depicted in Fig. 3. Different curves correspond to the different "spatial sizes" 
R and different "lifetimes" r of a hot pion system at the temperature T = 180 MeV. Notice that 
the production rate in a finite small pion system differs from the rate in a infinite pion gas (solid 
curve) where pion in-states can be taken as plane waves. The deviation bigger when the parameters 
R and r are smaller. Of course, this is a reflection of the uncertainty principle which is realized by 
the presence of the distribution \p{K — P)p as the integrand factor in 0. Basically, the presence 
of the form-factor of the multi-pion system will result in a broadening of the rate for small invariant 
masses M ^ 800 MeV/c^ which is wider at the smaller parameters R and r. This seems natural 
because the quantum fluctuations of the momentum are more pronounced in smaller systems. We 
emphasize as well that the behavior of the curves in Fig. 3, which correspond to a finite system, has 
a similar tendency to the CERES data Pfl |2] . 
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Fig. 3. Rates of electron-positron (left panel) and muon-muon (right panel) productions in pion-pion annihilation in 
a small finite system, T = f 80 MeV. 

For comparison, we present in Fig. 4 the results of evaluation of the rate dR^f+^_ / dMdy of electron- 
positron pair production in quark-antiquark annihilation in a QGP drop. The evaluations was made 
under the same assumptions as for pion-pion annihilation. As in the previous case, one can see 
an increase in the rate with a decrease in the invariant mass up to two electron masses. This 
real threshold is close to the total mass of annihilating quarks M = 2mq « 10 MeV/c^. The 
different behavior of the rate for small parameters R and r in the region of small invariant masses 
M ^ 500 MeV/c^, as compared to the rate for infinite parameters R = oo, r = oo, is due to an 
increase of quantum fluctuations which are evidently bigger for a smaller size of the QGP drop. 

One more important result is worth to note: the enhancement of the dilepton production rate 
for the low invariant mass region is much more sensitive to the variation in the spatial size of a 
many-particle (pion, quark) system than to that in the system lifetime (see Fig. 4, right panel). This 
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Fig. 4. Electron-positron production rates in quark-antiquark annihilation in a hot QGP drop, T = 180 MeV. 

fact can be explained with the use of the particle density taken in a particular form ((7j). Indeed, for 
the system which is confined in three spatial dimensions, each of the three spatial densities in (0), for 
instance exp {—x1/R^) with i = 1,2, 3, corresponds to one integration in ((T)) over a component of the 
pion-pair total momentum Ki. Hence, each of the three integrations on the r.h.s. of (jH]) involves a 
factor which is responsible for quantum fluctuations in its own spatial dimension; for instance, one of 
the three factors is (27r)^/^/2-exp {—R^Kf). Meanwhile, the spectral function which is responsible for 
the system lifetime finiteness, for instance (27r)^/^r ■ exp {—t'^Kq) as in our approach, emerges as an 
integrand factor just once during the integration over Kq. That is why the 3-dimensional variations 
of the spatial size of a many-particle system, as seen in Fig. 4, influence the shape of the dilepton 
production rate much stronger than the variations of the lifetime of a small system. 
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